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Abstract. Initially motivated by their relevance in foundations of quantum mechan¬ 
ics and more recently by their applications in different contexts of quantum informa¬ 
tion science, violations of Bell inequalities have been extensively studied during the last 
years. In particular, an important effort has been made in order to quantify such Bell 
violations. Probabilistic techniques have been heavily used in this context with two 
different purposes. First, to quantify how common the phenomenon of Bell violations 
is; and secondly, to find large Bell violations in order to better understand the possibil¬ 
ities and limitations of this phenomenon. However, the strong mathematical content 
of these results has discouraged some of the potentially interested readers. The aim of 
the present work is to review some of the recent results in this direction by focusing 
on the main ideas and removing most of the technical details, to make the previous 
study more accessible to a wide audience. 


Introduction 

Bell inequalities have attracted much attention in the last years. Their original inter¬ 
est as a key tool in the study of foundations of quantum mechanics has been nowadays 
surpassed by the relevance of these inequalities in different contexts such us quantum 
cryptography, communication complexity protocols and generation of trusted random 
numbers. In addition, Bell inequalities have been shown to be intimately related with 
some problems in computer science and operator algebras theory, capturing in this way 
the interest from those communities. Given the great importance of probabilistic tech¬ 
niques in those fields as well as in different areas of quantum information, it is not 
surprising that they are also very useful in the context of Bell inequalities. In fact, ran¬ 
dom constructions have been a key tool to solve some questions which had remained open 
for a long time in the field. However, despite their potential usefulness, these techniques 
and results are still far from being considered natural by many people working on quan¬ 
tum nonlocality. The aim of the present work is to review some of the most important 
results in the context of Bell inequalities for which probabilistic techniques have played 
a crucial role. Here, we will focus on the main ideas without paying attention to the 
technical details with the hope that this makes the previous works more appealing for 
the non-experts. Hence, this work must not be understood as a general survey on Bell 
inequalities, for which the reader can find excellent references in (TSj, |T6j . [46] . [45] . In 
particular, we will deliberately skip some standard topics such as connections with other 
areas, physical applications of the results and so on, with the upside of going directly to 
the important points of our discussion. 
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Let us start by introducing the basic objects of study. Bell inequalities were first con¬ 
sidered by Bell in [8] as a way to clarify an apparently metaphysical dispute on the com¬ 
pleteness of quantum mechanics as a model of Nature arising from the work [23] . Given 
two spatially separated quantum systems, controlled by Alice and Bob, respectively, and 
described by a bipartite quantum state p, Bell showed that certain probability distribu¬ 
tions obtained from an experiment in which Alice and Bob perform some measurements 
x and y in their corresponding systems with possible outputs a and b, respectively, 
cannot be explained by a classical modest More precisely, if P = (P(a, 6|x, y))^ A a fe=1 
denotes the probability distribution obtained in such an experiment, we say that P is 
a classical (or local) probability distribution if it can be written as 

P(a,b\x,y) = I P UJ (a\x)Q u (b\y)dP(u:) 

J Q 

for every x, y, a, b, where (12, P) is a probability space and for every oj € 12 we have 
P u (a\x) > 0 and Y.a P ^( a \ x ) = 1 f° r a i x ( res P- Qcj(b\y) > 0 and ^2 b Qu(b\y) = 1 
for all b,y). We denote the set of classical probability distributions by C. It is easy to 
see that C is a polytope (convex set with a finite number of extreme points) and the 
inequalities describing its facets are called Bell inequalities. On the other hand, we say 
that P is a quantum probability distribution if 

P(a, b\x, y) = tr(£% <g> Fyp) 

for every x, y, a, b, where p is a density operator acting on the tensor product of two 
Hilbert spaces Hi^H? and (E x ) Xja , ( Fy) y ,b are two sets of operators representing POVM 
measurements acting on H\ and H2 respectively. That is, El) > 0 and Y) a ^x = Iffi 
for all a, x (resp. Fy > 0 and F,j) = 11 # 2 for all b,y). We denote the set of quantum 
probability distributions by Q. It is easy to see that this set is convex and it verifies 
jC C Q. As it was shown by Bell, the converse inclusion fails; equivalently, there exist 
quantum probability distributions violating some Bell inequalities. 

In fact, in the beginning of this theory a slightly simpler scenario was considered. In 
the particular case where Alice’s and Bob’s measurements are binary (that is, a x ,b y = ±1 
for every x, y) one can consider the joint correlation: 

lx,%, = E[a x • by] = P( 1, l|x, y) + P(- 1, -l\x, y) - P(- 1, l\x, y) + P( 1, -l|x, y) 

for every x,y = 1 By plugging the definition of classical (resp. quantum) 

probability distribution in the previous expression one justifies the definition of classical 
(or local) correlation matrices as those which can be expressed in the form 

1x,y — / (cu)Hy((u)dP(cc) 

Jil 

for every x, y, where (12, P) is a probability space and for every cu € 12 we have A x (uj) € 
{—1,1} for every x (resp. B y (ui) € {—1,1} for every y); and quantum correlation matrices 


^Formally, Bell talked about a local hidden variable model (LHVM). 

^Note that P is not a probability distribution itself. For every x, y fixed we have that (P(a, b\x, y))a : b=i 
is a probability distribution. However, it is standard to use this terminology. 
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as those of the form 


T x,y — \ x ® Byp)., 

for every x, y , where p is a density operator acting on the tensor product of two Hilbert 
spaces H\®H 2 and ( A x ) x , ( B y ) y are self-adjoint operators acting on H\ and H 2 respec¬ 
tively verifying max Xi j / {||H a: ||, ||B y ||} < 1 . Since the different contexts will be clear along 
the work, we will also denote by £ and Q the set of classical and quantum correlation 
matrices respectively. Note that in order to consider this simpler context we made two 
simplifications: We restrict to binary measurements and we only considered the joint 
correlations of Alice’s and Bob’s measurements (and not the marginals). In particular, 
£ is again a polytope defined by its facets, now called correlation Bell inequalities , con¬ 
tained in the convex set Q. Bell’s work actually showed that this inclusion is strict even 
for the simplest case x,y = 1,2 (which implies the result for probability distributions). 


In this work it will be useful to understand Bell inequalities in a slightly more gen¬ 
eral sense as dual objects of probability distributions. In particular, any real tensor 
M = (Mx’,y)xy X a b = i (resp. real matrix M = (M X:y )^ y=1 ) defines a Bell inequal¬ 
ity (resp. correlation Bell inequality) by considering the dual action on probability 
distributions (resp. correlation matrices): ( M,P) = Yl x ya b Mx’,yP( a ,b\x } y) (resp. 
(M, 7 ) = YlxyM x ,ylx,y)- Given M, we will denote its classical value and its quan¬ 
tum value respectively by 

u(M) = sup{|(Af,Z)| : Zg£} and w*(Af) = sup {|(M, Z)\ : Z G Q). 


This picture allows us not only to describe Bell inequality violations: 


M, LV{M ) := 


w*(M) 

w(Af) 


but also to quantify these violations by means of the above quantity LV. As we will 
explain below, this quantification is crucial if one wants to understand the possibilities 
and limitations of Bell inequalities. The most famous correlation Bell inequality, the 
CHSH inequality [19], is given by the 2x2 matrix M defined by M 2 ; 2 = — 1 and 
M x y = 1 otherwise, for which one can check that cj(M) = 2 and w*(M) = 2\/2. Thus, 
LV(M) = y/2. 


Although the previous descriptions only consider the bipartite case, we can extend 
all the previous definitions to the multipartite setting straightforward. The multipartite 
scenario will be very important in our work because of two reasons. First of all, even 
in the tripartite case we will observe new phenomena which cannot be found in the 
scenario of two parties. Secondly, considering n parties introduces a new parameter in 
the problem and, as we will see, the answer for different questions can strongly depend 
on it. 


Once we know that violations of Bell inequalities exist, it is natural to wonder how 
common this phenomenon is. As we explained previously, there are two key objects 
in our picture, Bell inequalities (resp. correlation Bell inequalities) and probability 
distributions (resp. correlation matrices), which lead to different questions. On the one 
hand, one can study the probability of having LV(M) > 1 if we pick a Bell inequality 
(resp. correlation Bell inequality) M at random or, even more, which value LV(M ) we 
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should expect if we follow this procedure. On the other hand, one can pick a quantum 
probability distribution (resp. quantum correlation matrix) Q at random and ask how 
likely the event Q (f C is. However, two main obstacles appear when studying these 
questions. 

The first problem is that it is not so clear what picking these objects at random 
means. For that, one has to consider a probability measure on the corresponding sets 
and that is not always easy nor natural. In fact, the sampling procedure is particularly 
tricky for quantum probability distributions (resp. quantum correlation matrices) since 
here one has to impose a certain structure which is not required when one samples Bell 
inequalities. Regarding the definition of quantum probability distributions and quantum 
correlation matrices above one could think about sampling states and measurements 
independently. However, this procedure presents several problems such as deciding the 
dimension of the corresponding Hilbert spaces, the lack of a natural way to sample 
general measurements and so on. Some works have also considered the problem of fixing 
one of these objects (either the state or the measurements) and sampling the other one 
at random. In fact, considering random measurement seems to be an interesting problem 
from the experimental point of view since it is directly connected to the absence of a 
shared reference frame in Bell experiments. 

A second problem is that the quantities u(M), u ;*(M) and LV(M ) are in general 
very difficult to compute. The study [J6j performed by Tsirelson clarified the context of 
bipartite correlation Bell inequalities. In particular, he gave an alternative description for 
the quantity u*(M) which allows to understand the quantity LV 2 (N) = sup {LV(M) : 
M = (M XiV )£ ija via the so called Grothendieck constant. However, the situation is 
more intricate in the case of general bipartite Bell inequalities as well as in the tripartite 
correlation case. If we define LV 2 (N, K ) = sup {LV(M) : M = (M“’^)^ y A a fe=1 }, it is not 
clear how large this quantity can be as a function of N and K and the same happens for 
the quantity LV^{N) in the tripartite correlation scenario. Hence, some questions should 
be answered before trying to understand the probabilistic behavior of Bell violations in 
these contexts. Interestingly, probabilistic techniques are also very important for this 
purpose since they have been shown to be a useful tool to study how large the quantities 
LVs(N) and LV^iV, K) can be. On the other hand, beyond its importance from the 
previous point of view, looking for large Bell violations is an interesting question itself, 
since the quantity LV can be understood as a quantification of how better one can solve 
certain tasks if quantum resources are used instead of classical ones (see [26], [28], [36] 
for more information). 

Finally, let us remark that, despite the restriction imposed in the title of this survey, 
the number of works on Bell inequalities using probabilistic tools is very large. Here, we 
will focus on those results studying asymptotic behaviors. That is, those results describ¬ 
ing the probabilistic nature of a problem when (some of) the parameters go to infinity. 
With this in mind, let us remind the reader the standard asymptotic notation, which 
will be constantly used in this work. Given two nonnegative functions /(n) and g(n ) on 
natural numbers, we will write f(n) = 0(yg(n)) (resp. /(n) = fl(g(n))) if there exist a 


■Hhe subscript 2 in LVi{N) denotes that it is in the bipartite case. 
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constant C and a natural number uq such that f(n) < Cg(n ) (resp. f(n) > Cg(n)) for 
every n > no- On the other hand, we will write f(n) = o(g(n)) if hnu^oo f(n)/g(n ) = 0. 

The survey is organized as follows. Section [T] is devoted to explaining some results 
about the probabilistic nature of quantum nonlocality in the bipartite correlation case. 
As we will see, the deep understanding of this situation thanks to Tsirelson’s work allows 
to study the problem in a very natural way. In Section [2] some results about lower an 
upper bounds for tripartite correlation Bell inequality violations will be reviewed. This 
will lead us to the first examples of unbounded violations of tripartite correlation Bell 
inequalities. Section [3] will deal with some results about the probabilistic nature of mul¬ 
tipartite Bell inequalities in some particular cases as a function of the number of parties 
n. Finally, in Section [I] some recent results about general bipartite Bell inequalities and 
the quantity LV 2 (N,K) will be reviewed. In particular, we will explain some random 
constructions in this setting and compare them with some new results from computer 
science. 


1. Bipartite correlation Bell inequalities 

As we mentioned in the Introduction Bell’s work showed the strict inclusion C £ Q (see 
also m)- This picture was completed by Tsirelson in the case of correlation matrices, 
by showing that the set Q is not much bigger than the set C. More precisely, one has 
£ ^ Q ^ where 1.67696... < Kq < 1.78221... is the real Grothendieck constant^. 

The following result is the standard statement of the corresponding theorem. 

Theorem 1.1 (Tsirelson). 

(1.1) LV 2 (N) := su p{ ^^ : M = ( M X ,y)x, y =i } < k g for every N. 

In fact, the Grothedieck constant can be defined by Kq := sup N LV 2 (N). Theorem 
o is a consequence of Grothendieck’s inequality and a result proved by Tsirelson [36j 
which states that 7 = ( 7 x,y)xy=i is a quantum correlation matrix if and only if there 
exist a real Hilbert space H and unit vectors iq, ■ ■ ■ ,un,v 1 , ■ ■ ■ , ujv in H such that 

(1.2) 7 x, y = {u x ,Vy) for every x, y = 1, • • • , TV. 

In [SB ] the author posed the open question of whether a similar result to Theorem 11.11 
holds in the tripartite case. This is related to the lack of Grothendieck’s inequality for 
trilinear forms and we will go over this question in the following section. 

In [T] the authors tackled the question of how likely it is for a random bipartite 
correlation Bell inequality M to verify that the quotient between and uj{M) is 

strictly larger than one. In order to study this problem, one first needs to define a way 
of sampling these inequalities. In [T], the authors considered random N x N matrices 
M sampled from {—1,1}^ with respect to the uniform measure. In fact, although the 
authors focused on sign matrices, the same techniques can be applied to study more 

4 The exact value of the Grothendieck constant is still unknown in both the real and the complex case 
(see m for the most recent progress). 

^The work [T] deals with bipartite XOR games, but the problem is completely equivalent. 
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general random matrices, such as real gaussian matrices. The main result in p] states 
as follows. 

Theorem 1.2. If M = j is a random matrix sampled from {—1, l }^ 2 with 

respect to the uniform measure, then 

lim p(m : 1.5638... > U > 1.2011...) = 1. 

N-> oo l w(M) J 

Theorem 1 1.21 implies that for almost any correlation Bell inequality its quantum value 
is strictly large than its classical one. On the other hand, according to Theorem 11.11 and 
the comments below, the previous equation also tells us that with probability tending 
to one these inequalities are not so close to the optimal value Kq € [1.67,1.79] for the 
quotient between uj*(M ) and uj(M). Let us also mention that, although very natural, the 
way of sampling bipartite correlation Bell inequalities considered in [T] is very particular. 
Indeed, this procedure does not consider all bipartite correlation Bell inequalities since, 
for instance, all the entries of the matrices have the same absolute value 1. One can 
think about some other ways of sampling them by considering for instance a bipartite 
correlation Bell inequality M £ P^ as an element in the unit sphere with respect to the 
norm ||M|| = Yl x y =l \M x ,y\- Then, one can naturally define a probability distribution on 
this sphere. Preliminary calculations suggest that similar results to Theorem 11.21 could 
be obtained in that context. 

Since the most important part of Theorem 11.21 is the lower bound for the quotient, 
we will just briefly explain how to prove that estimate. To this end, the authors first 
show that liniAr-^oo P{M : uj(M) < (1.6651... + o(l))./V 2 j = 1. This is obtained as an 
application of the Chernoff bound to the random variable Ylxy =l e x,yt x s y for a fixed 
choice of signs t x ,s y = ±1, x,y = l,--- ,1V; and a counting argument to consider the 
2 2N possible choices of signs. In order to prove the estimate liniAr-^oo P{M : cj*(M) > 
(2 — o(l))7V5} = 1, from where one obtains the lower bound in Theorem 11.21 the 
authors used a clever construction based on the Marcenko-Pastur law, which describes 
the behavior of the singular values of the random matrix M. If we call L and R the 
n x m matrices whose columns are respectively the left and right singular vectors of the 
matrix M associated to its m largest singular values, basic linear algebra shows that 

N m 

(1.3) ^ [ M x ,y('U / xi Vy) — ^ ^ 

x,y= 1 i= 1 

where the m-dimensional vectors u x s and v y s are the rows of L and R respectively and 
(Aj)’Ti are the corresponding singular values. On the other hand, the Marcenko-Pastur 

law [32] tells us that if one defines the function f{s) = ^ yj^ — 1 dx on [0,2], for 

every e > 0 the number m of singular values verifying Aj > (2 — e)y/N belongs to the 
interval [(/(2 — e) — o(l)) N, (/(2 — e) + o(l)) IV] with probability tending to 1 as N goes 

to infinity. Hence, the quantity (11.31) is lower bounded by (2 — e)(/(2 — e) — o(l)) IV 2 with 
probability tending to one. The technical part of the proof in |Tj consists of adapting 
the Marcenko-Pastur law to show that one can assume that for a given 5 > 0, all the 
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previous vectors u x s and v y s have norm lower than or equal to \J /(2 — e) + 5 with 
probability tending to onc0. Therefore, by considering the normalized version of the 
previous vectors u x , v y , one obtains that the quantum correlation 7 = ((n x , ^ y )) xy=1 
verifies 


N 


U 


( M ) > ^ M x ,ylx,y ~ 


x,y=l 


f(2-e) + 5 


E A < 

i— 1 


> 


(2 — f ) (/(2 — c) — o(l)) , r> 

/{2-e) + J 


Since e and 5 can be made arbitrarily small one obtains the desired estimate. 

3 

Interestingly, the authors show in [T] that (2 + o(l))JV 2 is also an upper bound for the 
quantum value so it is optimal. However, the best lower bound for the classical 

value to(M) obtained in [T] is far from the upper bound explained above (see [T] Section 
4] for details). Obtaining the exact asymptotic value of uj(M) looks a difficult open 
problem, which seems more interesting from a mathematical and computer scientific 
point of view than from a physical perspective. 


In | 23 ] the same problem was considered from the correlation matrices point of view. 
That is, if one picks a quantum correlation matrix at random, which is the probability 
that it is nonlocal? The equivalent reformulation (11.21) of a quantum correlation gives 
a natural sampling procedure which avoids most of the problems that we mentioned 
in the Introduction: one can pick the vectors u\, ■ ■ ■ , un, v\, ■ ■ ■ ,vn independently and 
uniformly distributed on the unit sphere of IR m . It is well known that this is exactly 
the same as sampling independent normalized m-dimensional real gaussian vectors. It is 
very easy to see that if one fixes any finite m , the probability that a quantum correlation 
matrix sampled according to the previous procedure is nonlocal tends to one as N tends 
to infinity (see [241 Section 2] for details). However, this kind of sampling does not say 
too much since the set of quantum correlation matrices of order N which can be obtained 
with a fixed m is very small. The interesting case is that where m and N are of the 
same order. The main result in [241J gives indeed and answer to the considered problem 
as a function of a = S: 


Theorem 1.3. Let u\, ■ ■ ■ ,un,v\, ■■■ ,vn be 2 N vectors sampled independently accord¬ 
ing to the uniform measure on the unit sphere of R m and denote by 7 = ((ui,Vj))^j =1 
the corresponding quantum correlation matrix. Then, if we denote a = ^ we have 

a) If a < ao ~ 0.004, then 7 is nonlocal with probability tending to one as N tends to 
infinity. 

b) If a > 2, then 7 is local with probability tending to one as N tends to infinity. 


There is a considerable gap between ao and 2. In fact, one should not expect this result 
to be optimal. This is because the proof of part a) above is based on an approximation of 
gaussian matrices by orthogonal ones and a use of the main result in [T], where N must 
be artificially larger than m. However, the important point of the previous statement is 
that it shows a nontrivial phase transition for the nonlocal properties of 7 as a function 

®In fact, the modified Marcenko-Pastur law proved in [T] Theorem 3] gives the probability for each of 
these vectors to have norm larger than y/f (2 — e) + <5. Then, the authors proved that Eq. (11.31) is not 
affected if one rules out those vectors with large norm. 
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of the parameter a = Clarifying the case a = 1, studying any possible relation 
between a and Kq and reducing the gap appearing in Theorem 11.31 are proposed in [23] 
as open questions. 

In order to understand how part a) above can be proved, let us go back to Theorem ll.2l 
and to reformulate it in a different way. Let us assume that we sample two independent 
N x N orthogonal matrices U and V according to the Haar measure. Then, if we 
call Ui (resp. Vj ) the normalized vector formed by the first m coordinates of the i- 
th row of U (resp. j-th row of V), then the quantum correlation 7 = ((itj, 
is not local with probability tending to 1 as JV goes to infinity. Here, m is the one 
considered in Eq. m above. Indeed, this can be deduced from the fact that if we 
consider the singular value decomposition of a real random gaussian matrix G = U DV 
(for which Theorem 11.21 can be stated exactly in the same way), the matrix U and V are 
independent Haar distributed orthogonal matrices (see [24l Propostion 2.1] for details). 
However, this fact cannot be used directly in [ 23 ], since in the problem considered there 
one must start sampling gaussian matrices G\ and G 2 instead of orthogonal matrices 
U and V. More precisely, according to the comments above, the problem considered in 
[ 23 ] consists of studying the nonlocal properties of the correlation ((g*, hj))^. =v where 
gi (resp. hj ) are the normalized vectors formed by the m first coordinates of the i-th 
row of an N x N gaussian matrix G (resp. j -th row of an N x IV gaussian matrix H ). 
However, part a) above can be obtained by invoking [25) Theorem 1.1], which shows that 
there is a coupling between real gaussian matrices G and Haar distributed orthogonal 
matrices O such that the norm sup i=1 N \\f™(g - Vivo)|| is controlled, where here 
FJ n {G — VNO) is the i-th row of the matrix G — '/NO truncated to its first m entries. 
Indeed, a suitable control on the previous norm allows the authors in [24] to replace 
the nonlocal correlation 7 introduced above from two orthogonal matrices by another 
nonlocal correlation 7 = ((gi, hj)). ._ v where the gi s and h/s are now m-dimensional 
normalized real gaussian vectors (see [ 231 Theorem 2.3] for details). 

Part b) of Theorem 11.31 can be obtained by using classical Banach space techniques. 


2. Tripartite correlation Bell inequalities: Unbounded violations 

2.1. Extensions of Tsirelson’s result to the multipartite setting. In order to 
study Bell violations in the multipartite case and regarding the importance of the GHZ 
state in the setting of two parties, it seems very reasonable to consider the generalized n- 
partite d-dimensional GHZ state: \ip) = A= /2i=] On the other hand, any possible 

strategy looking for unbounded Bell violations should definitely exploit the absence of 
the Grothendieck inequality in the multilinear framework. However, one should precise 
this statement a little bit. Grothendieck’s inequality can be stated in many equivalent 
ways (see Ell Page 172]) and it turns out that, in the multilinear case, the corresponding 
generalizations of those statements are not equivalent anymore. In fact, although several 
of the possible extensions of Grothendieck’s inequality to the multilinear setting have 
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been proved to be false, a few of them remain valid in the new context (see [9], [TO] . [18], 
[33]). The following generalization of Grothendieck’s inequality was proved in [9], @30 


Theorem 2.1. Let n,N > 2 and d be positive integers, let T = (Tj 1; ... ,j n )^ ... i)]=1 be 
a real tensor and xq, • • • ,Xi n elements in the unit ball of a complex Hilbert space H of 
dimension d for every i\, ■ ■ ■ , i n = 1, • ■ ■ , N. Then, 

N 


E 


L U»— 3 In V^l 5 


(Xi 


,Xi 


3n 5 /-i . 

<2—K c g -u(T), 


'i'li’" fin — 1 

where (xq,--- ,Xj n ) = Ylt =i x h (&)''' x i n (&) and Kq is the complex Grothendieck con¬ 


stant. 


In [SB; Theorem 11 ] the authors used another version of Theorem 12 .II to show that the 
largest Bell violation achievable by three partied sharing a d-dimensional GHZ state is 
upper bounded by 4\/2 Kq, independently of the number of inputs N and the dimension 
d. This was the first result proving that a nontrivial family of states can not give 
large Bell violations and it can be seen as a generalization of Theorem 11.11 In [36] the 
authors posed the problem of whether an analogous result could be proved for Schmidt 
states: \i/j a ) = Yli=i a i\i)® n , where a = (cu)f =1 verifies |a*| = 1 ; and provided 

in addition a direct connection between such a problem and an open question in the 
context of operator algebras. Two years latter this question was answered in P31 by 
using a surprisingly easy argument. The authors in m realized that the upper bound 
for the GHZ state can be obtained in a straightforward manner from Theorem 12.11 and 
extended the result to Schmidt states by using a nice expansion of each of these states 
in terms of non-normalized GHZ states (see m Theorem 1]). 

Theorem 2.2. Let T be an n-partite correlation Bell inequality. Then, for every n- 
partite quantum correlation 7 constructed with the state \if a ) = Ylo-i we have 

that {T, 7 ) < 2 2 K g ■ oj{T), independently of the number of inputs N and the local 
dimension d. 

In [13] an exhaustive study of Carne’s extension of the Grothendieck inequality to the 
multilinear case m was performed to conclude that an analogous result to Theorem 12.21 
can also be stated when the parties share a clique-wise entangled state (see m Theorem 
2]). This implies in particular that for tripartite correlation Bell inequalities, the amount 
of Bell violation achievable by an arbitrary stabilizer states is uniformly bounded. 


2.2. Unbounded violations of tripartite Bell inequalities. The previous results 
rule out most of the candidates one would first study in order to find large Bell violations 
in the multipartite setting. Then, as in many other contexts, it becomes natural to study 
the behavior of random states. A standard way of sampling random pure quantum states 
1^) € C d is by using the uniform measure on the unit sphere of the Hilbert space 

’-7 

In fact, we state here a slightly modified version for T real, which involves a modification in the constant 
(see 1131 Theorem 9] for details). 

®The result can be generalized to n parties straightforward. 
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or, equivalentl)@, using gaussian random variables | ij>) = Yli=i 9iV) G C d . When we 
are dealing with more than one system this last notation is usually extended to | if) = 
JT ... t =i 9ii,-,in\h ' "in) € (C d )®”. In the case of two and three systems it is also very 
common to use random unitaries sampled according to the Haar measure in the unitary 
grupo V d . We write \if) = 1 U(i,j)\ij) € (C d )® 2 or \if) = J2i,j,k=i u iti> k )\ijk) G 

(<C d )® 3 , where in the last expression the unitaries (Ui)f =1 are sampled in the cartesian 
product Urf. The unitary approach has been extensively used in quantum channel 
theory and entanglement theory. 

In [36] it was proved for the first time that, in contrast to the bipartite case, there are 
tripartite correlation Bell inequalities which lead to unbounded violations. 

Theorem 2.3. For every dimension d £ N, there exist D € N, a pure state \if) € 

rl n D / \2 rf2 2 D2 2 d2 

C a ® ® and a Bell inequality T = (T ^*.) i ■ ^ =1 ’ such that the violation by \'ip) 

on such an inequality is 

This theorem implies that there does not exist a uniform constant C such that 
LVs(N) < C independently of the number of inputs N ; giving in this way a nega¬ 
tive answer to the question posed by Tsirelson in [36j. It is advisable to extend the 

previous definition to LV^N^d) = sup {: T = (T x , y ,z)x, y ,z=i}i where here w* d (T) 
denotes the quantum value of T when we only consider tripartite quantum correlations 
constructed with quantum state of local dimension d. Then, it is important to point 
out that, in order to obtain unbounded violations, one must allow to increase both the 
number of inputs N and the dimension of the Hilbert space of the system d. Indeed, 
if we fix one of these parameters, then the amount of violation is upper bounded by a 
constant depending on it. 

Theorem 2.4. The following upper bound holds: 

LVz(N , d) = O (Vk) , where k = min | IV, d}. 

The upper bound as a function of d was first proved in j36] (see m Theorem 3] for 
an alternative proof). In fact, one can state a stronger result since it suffices that only 
one party has local dimension d. This tells us that Theorem 12.31 is optimal in the local 
dimension d. The (easier) upper bound as a function of N can be found in m Theorem 
2] and it also admits an extension requiring that only one party has N inputs. In |14l 
the authors generalized the previous theorem to n parties by providing the upper bound 

n—2 

0(d ^~) whenever the share state of the parties is restricted to have local dimension d 

n—2 

on at least n — 2 players and to 0(N ^~) whenever the inputs in the correlation Bell 
inequality is at most N for at least n — 2 players. 

The key point to prove Theorem 12.31 is the use of random states, showing once more 
that these states exhibit unexpected extremal properties. The proof of the theorem relies 
on hard techniques from operator spaces and it is highly nonconstructive. In particular, it 
does not provide neither an explicit (nor even probabilistic) form of the Bell inequality T 

®One needs to normalize the gaussian state to have the same distribution. 
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attaining such a violation nor any control on the dimension D appearing in the statement 
of the theorem. Since we will explain below a more recent result improving Theorem 
m we will not say too much about the proof of the previous result. Instead, we will 
present some basic ideas in order to be able to discuss the analogies and differences with 
the later proof. 

The first idea in the proof of Theorem 12.31 is to reduce the problem to work with 
Hilbert spaces. Indeed, the initial problem consists of comparing two difficult to handle 
quantities w(T) and c v*(T) for a given tensor T = {Ti,j,k)fjk= i- Let us consider another 
tensor S = (•SVj.fc)”) fc=1 as an element in C m < 8 > C m C m and define two new quantities 
(norms): 


\\s\U 


jrn = SUP 


^ ^ bj C-k 

ij y k= 1 


: |iMr.iii2.iife)?Liii 2 .iife)r„ii2<i} 


and 


m 

II'S'IU = sup | I ^2 Si,j,kAi®Bj®Ck ■■ ||(A)™i||rc, \\(Bj) 1 jLi\\RC, \\{Ck)T=i\\RC < l}, 
i,j,k =i Md3 


where for a sequence of complex numbers (zj) 7 JL 1 , || (jstj) y!L 1 II 2 is the euclidean norm and 
for a sequence of matrices (Zj)™ =x c Md we define 


(2.1) ||(Z J )”l 1 || i?C7 = max{||^Z J zj||^||^zjz j ^}. 

i= 1 i =1 

These quantities can be understood as a hilbertian version of the values w(-) and <*;*(•). 
On the other hand, for every tensor S one can construct another tensor T = {^Ti,j,k)fj k=i 
with N = 2 m for which, up to a universal (known) constant, 

u*(T) __ ||S||» 

U} ( T ) \\S\Uz i ® c ez i ® e ez i 

Here, we should point out that one could explicitly construct both the Bell inequality 
T and the observables to be used to compute ui*{T) from the elements S , (H*)™ j, 
( Bj )jli, (C’fc)fcLi used to compute ||S||*. In order to find a tensor S for which the 
quotient between ||5||* and is large, one considers the random state 

IV’) = m ^Tj k=i Ui(j, k)\ijk) € (C m )® 3 . It follows from well known results in random 
matrices that ('01V’) — 1- Then, by doubling indices one can naturally consider the 
element 5 in C™ 2 <8 C™ 2 <g> C™ 2 defined by = {k\U} r \j){k'\ul\j') and the 

matrices A iy = -A=\i){i'\, B jtj > = C k ,k' = ^ \k)(k'\ for every i, k, k’ = 

1, ■ ■ ■ ,m. It is not difficult to check that for these matrices the quantity (12.11) is lower 
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than or equal to one, which implies 
1 


(2.2) ||S||* > — U) ^2 S it i>jj' tk ,k'\i)(i' | ( 8 ) \j){j'\ (8 \k)(k'\ i/>\ 

7712 i,i',j,j',k,k '=1 

1 m - 

= — J2 tr ((pi®ul){ul r ®Uif) \ = — tr(l Mm2 ) = 

m 9 . ' ' jji 2 




According to the previous comments, we would finish the proof if we could show that 
||iS|L m 2 „ m 2 2 < 1. Unfortunately, it turns out that this estimate is not true. The 

key result [36l Proposition 20] shows that for every e > 0 and for every m € N , there 
exist N(e,m) € N, some unitary matrices Ui € Mn and some matrices H^y € 2 for 

every i, i! = 1, • • • , m such that if we define Fiy = C/* r <8 t/t + H,y and Siy-jj/yy = 
(kk'\Fiy\jf) , we have that 


(2.3) 

ll^ll^ 2 ® e ^ 2 ®^f 2 ~ an< ^ 


-j-ptr(Hiy(Ui®Ui'))\ < e for every i,i' = !,-■■ ,m. 


This means that, while the value of the norm II • ||. m 2 „ „ N 2 is smaller for the new 

inequality S , the modification with respect to the previous inequality does not affect 
essentially to the estimate ( 12 . 21 ) . Hence, the result follows. As the reader can guess, it 
is precisely in the proof of the existence of these highly nontrivial matrices H^y where 
the explicitness of the result is lost. 

In [37], [38j the author refined the preceding proof by using random gaussian matri¬ 
ces instead of random unitaries. The problem is again reduced to separate the norms 
and US’!!* but the technical proof to find the previous unitaries is replaced 
by the use of a previous known result (EH Theorem 16.6]), which seems to be due to 
Steen Thorbjprnsen. This result guarantees the existence of a family of N x N gaussian 
matrices for a certain N, for which the construction explained above can be 

done directly. The main advantage of this approach is that one can use directly these 
gaussian matrices and the matrices H,y are not needed anymore. Actually, one can fol¬ 
low the estimates behind these results and replace the previous transformation from the 
tensor S to T by a more sophisticated one (allowing N ~ n 2 ) to obtain a Bell inequality 
T = (Tj jy ) - j , k=1 I which can give violations of order y/N by using a quantum state in 
(8 C ^ 2 (8 C ^ 2 (see [Ml Remark 3.3] for details). 


2.3. Briet and Vidick’s construction. In the work [14] the authors gave another 
proof of Theorem 12.31 which considerably improved both the estimates on the parameters 
and the construction. 

Theorem 2.5. Let us assume that N = 2 3 for some natural number j. There exist a 
quantum pure states \ip) in C N <8 (8 C N and a Bell inequality T = (7ij,fc)fj 

such that the violation by \ip) on such an inequality is fl(\/lVTog - ' 5 / 2 IV). Moreover, the 
observables used in each party are tensor products of Pauli matrices. 
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More generally, for every N that is a power of 2 there exist an n-party Bell inequality 
T with N 2 inputs in each party and a state \ip) € (C Ar )® n such that they can lead 
to a violation £l((N log ~ 5 , where the observables used in each party are tensor 

products of Pauli matrices. 

According to Theorem 12.41 and the comments below it, the previous result is opti¬ 
mal, up to a logarithmic f'actoiPl. in the dimension of the Hilbert spaces and it is only 
quadratically off from the best upper bound as a function of the number of inputs N. 
The proof of Theorem 12.51 is probabilistic, being based again on the construction of a 
random Bell inequality which interacts properly with a random pure state via some suit¬ 
able observables. In fact, since the elements involved in Theorem 1231 are very simple, 
we will first explain them and a brief explanation at a more mathematical level will be 
discussed later. 

Let us consider the random state 

1 N 

( 2 - 4 ) l^) = jni| 9i,j,k\i,j,k), 

where g is the corresponding non-normalized state. A key point in m is the use of Pauli 
matricetO to show that if one considers a six indices tensor j y k fc'=D one 

can easily define a Bell inequality T = (Tp^q^p^q^r, where the inputs are indexed in 
the tensor product of j Pauli matrices and such that if the three parties use respectively 
the observables P , Q and R associated to the inputs P, Q and R, then it turns out that 

(2.5) ]T T PtQiR (<p\P < 8 >Q® R\<p) = N 3 {<p\S\tp). 

P,Q,R 

Here 5 = Si,i'-,j,j';k,k'\ijk)(i'j'k'\ is regarded as an element in M N 3. Indeed, 

this follows easily from the fact that the set Vj = {Pauli matrices}® 3 ' forms an orthogonal 
basis of Mjy with respect to the inner product defined as {A, B) = tr(AB^). Moreover, 
it is trivial to check that (P, Q) = N5p t Q for every P,Q& Vj. Then, by considering the 
set Vj < 8 > Vj <S) Vj one obtains an orthogonal basis of M N 3 , and for every element S in 
this space one can write (its Fourier expansion) 

5 = ^ ( S,P®Q®R)P®Q®R. 

P,Q,ReVj 

Hence, if one defines T = (Tp,q,_r)p,q,rg'P : , such that 

N 

Tp,Q,R — {S , P <55 Q ® R) — } " Siy-jji;k t k'Pi,i'Qj,j'Rk,k'i 

i,i';j,j';k,k'=1 

one has the desired property (12.51) . Note that in order to obtain this property one needs 
to double indices as it was made in the proof of Theorem 12.31 

^Pisier has shown in [38] that such a factor can be reduced to log - 2 N. 

^The only property used by the authors is that the tensor products of j Pauli matrices form an orthog¬ 
onal basis of Mjv formed by observables. Any other such a system would equally work in the proof. 
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On the other hand, let us assume that : Rj {—1,1} are optimal func¬ 
tions to compute w(T). By defining the hermitian matrices X = x(H)-P, Y = 

Sqs'Pj u (Q)Q-> Z = Y^ReVj C(R)R i n M n one can write 

u(T)= T PtQtWC {P)v(Q)aR) = {S,X®Y®Z). 

P,Q,R&Vj 

Then, the fact ||V ||2 = ||V ||2 = \\Z \\2 = N 3//2 implies that 

(2.6) uj(T)<N 9/2 sup ( S,X®Y®Z ), 

X,Y,ZeB(Hcrm(N)) 

where _B(Herm(A r )) denotes the unit ball of the hermitian N x N matrices with respect 
to the Frobenius norm. 

By looking at Eq. (ESI), a potential good choice for the tensor S is Siy-jj'-ky = 
9i,j,k9i’,j',k’ for every i,i' , j, f ,k,k'. Indeed, if one considers this element it is straight¬ 
forward to check that 

(2.7) oj*(T) > IV 3 (</?|S'|(^) = IVj|g||| ~ IV 6 with high probability over g, 

where the last estimate follows from well known results on gaussian variables. Therefore, 
the statement would follow if one proved that the quantity sup.Y,y,zeB(Herm(JV)) (^\ X <S> 
Y <S> Z) is upper bounded by N up to, maybe, some logarithmic terms. Unfortunately, 
it can be deduced from known results that the previous amount is 0(IV\/lV). The key 
point in the construction by Briet and Vidick is to remove some of the indices of S. 
More precisely, they consider the tensor S defined by = gij t k9i',j',k' whenever 

i / i', j 7 ^ j 1 , k / k' and Siy-jj'-k^k 1 = 0 otherwise. It is very easy to see that the 
corresponding values has the same order IV 6 , since the number of removed terms 

is negligible in the estimate (12.71) . On the other hand, the main result in m shows that 
the classical value of the new inequality does decrease by the previous modification: 

(2.8) sup (S,X (g> Z) < N log 5//2 N with high probability over g. 

X,Y,Z£B(Herm(N)) 

This estimate allows to obtain Theorem [23J Note that here, as in the proof of Theorem 
12.31 the definition of the Bell inequality T can lead to complex coefficients Tpq h = 
^Ai'-jj'-k k'=i r Yi,i'-,j,j'-,k,k'Pi,i'Qj,j'Rk,k'- However, it is trivial that either the real part or 
the imaginary part must lead to a similar violation up to a constant 2 . 

The proof of Eq. (12.81) is based on a technical e-net construction using a decomposition 
of elements in B(Herm(IV)) as linear combinations of normalized projections. In their 
proof the authors show that the corresponding estimate holds for such projections with a 
sufficiently good concentration so that, by applying a counting argument on the elements 
of the net, they obtain the result. 

A more careful look at the previous argument allows to see some analogies with the 

proof of Theorem l2.3l land the subsequent improvement in [38]). Indeed, one can see that 

2 2 2 

the proof by Briet and Vidick is also reduced to separate two norms in C m <8> <C m ~ <8> C m , 
being the first one again the e norm ||>S'||^ m 2 (gi prfi- However, in this case the second 
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norm is: 

m 

II-S1I3.3 = || ^ ^ ji k)(i ,j >^)||j^ v3 - 

The reader should note that, while in the proof of Theorem 12.31 the comparison between 
the norms consisted of computing the norm of the identity map on <8)f =1 C m when this 
space is endowed with the e and the || • ||* norm respectively, in the proof of Theorem 12.51 
one must compute the norm of the rearrangement map R : (g)f =1 C m —>• (8>^ =1 C m defined 
by R(\ii')\jj')\kk')) = \ijk)\i'j'k'). Note that doubling indices is an essential point to 
define this map. 

There are two key points to understand the improvement in Theorem 12.51 with respect 

to the previous results. The first one is that the transformation from the tensor S in 
2 2 2 

C m <S> C m (g) C m into the Bell inequality T = T(S) does not imply an increase in the 

number of inputs. That is, one has m = N. The second point is that the authors gave 

2 2 2 

an essentially optimal separation from the norms they considered in C m (g> C m (g> C m 
(see [38 . Section 2] for details). Finally, we should mention that Theorem 12.51 is not 
just an existence result, but the proof shows that such an estimate happens with high 
probability in the choice of the gaussian variables g. 

Theorem 12.51 can be understood as a result proving that tripartite correlation Bell 
inequalities give large violations with high probability when they are properly sampled. 
However, this sampling is rather artificial since, in particular, it implies doubling indices 
in the coefficients. Regarding SectionQ]in this survey, a next step in the problem could be 
to study the behavior of T = (Gi,.j,fc)fj fc=1 , where it is sampled from {—1, l}^ 3 according 
to the uniform measure. In fact, this is equivalent to consider T = fc=1 , a family 

of independent real gaussian variables. It is not difficult to see that E[w(T)] < N 2 
in this case. Any nontrivial estimate on E[w*(T)] would definitely be an interesting 
result. On the other hand, despite the sharpness of Theorem 12.51 there is still a gap 
with respect to the best known upper bound for the multipartite Bell violation as a 
function of the number of inputs. It would be interesting to know if one can indeed 
attain violation of order yfN by using a tripartite correlation Bell inequality with N 
inputs per party. Finally, any explicit (non-probabilistic) construction of a tripartite 
correlation Bell inequality leading to unbounded violations would be also very welcome 
by the community. 

3. Correlation Bell inequalities for a large number of parties 

3.1. Symmetric XOR games. In the works m, 0 and [lj Ambainis and coauthors 
studied the scenario of binary inputs correlation Bell inequalities with many parties 
T = (T Xl ,— , Xn )xi,- ,x n e {o,i}- I n this case, it is well known that the quotient between 
oj*(T) and oj(T) can be equal to 2 2 for some particular inequalities T (0, |33j). More¬ 
over, it is also known that such a violation is optimal in this setting [30], [50]. The 
reader should immediatly note that this problem is completely different from the one 
considered in the previous section, where unbounded Bell violations were proved by 
considering only a fixed number of parties n (say three) and increasing the number of 
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inputs N. In order to clarify the connection between the works Ei, 0 , m and the 
explanation in this survey, let us recall the reader that any correlation Bell inequal¬ 
ity T = {T Xlt ... tXn ) Xlt ... tXn verifying J2 xi ,-,x n \ T xi,-,x n \ = 1 can be trivially written as 
T = ( 7 r(aq, • • ■ ,x n )f(x i, ■ ■ ■ , ®n)) x x > where 7 r is a probability distribution over the 
set of inputs and / : Zj x ■ x Jv ^ {—1,1} is a function. This gives a correspondence 
between the classical (resp. quantum) value of a correlation Bell inequality and the 
classical (resp. quantum) bias (probability of winning minus probability of loosing) of 
an XOR game. The corresponding game is defined by the distribution n over the set of 
questions (inputs) and the predicate function, given by V{x\, ■ ■ ■ ,x n ,a\, ■ ■ ■ ,a n ) = 1 if 
and only if ai © • ■ • © a n = ^(/(aq, • • • ,x n ) + 1). Then, a symmetric XOR game is a 
correlation Bell inequality where n is the uniform probability distribution over the set 
of inputs and the function / is invariant under permutations of the parties^. Note that 
for the case of binary inputs aq = 0,1, a symmetric XOR game T can be described by 
means of a sequence of n+1 bits (Mo, Mi, • • • , M n ), where Mj = /(aq, • • • , x n ) whenever 
Y^i=i x i = j- Then, one can define a probability distribution on the set of these Bell 
inequalities by picking the (n + l)-bit string (Mo, Mi, • • • ,M n ) uniformly at random. 
The main results in |2j and 0 imply that for every e > 0 there exist positive constants 
01 (e), 02 (e) such that for a high enough n we have 


(3.1) 


Pj T : 


cO(T) 


€ 


Oi (e) \/log re, O 2 (e) ^J\og 



> 1 - e. 


Therefore, large Bell violations occur with probability very close to one on the set of all 
symmetric XOR games. On the other hand, since the violation 2 2 mentioned before is 
attained on symmetric XOR games, one sees that the violation ylogn is far from being 
optimal. 

The key point to prove Eq. m is a simplification of both quantities cj(T) and 
when T is a symmetric XOR games 0. On the one hand, in 0 the authors 
proved that in order to study uj(T) in this case, it suffices to look at n + 1 deterministic 
correlations 70, • • • , € C. These correlations are very easy to describe in terms of the 

strategy followed by the players to play the XOR game corresponding to T. Indeed, 
7^. is the correlation obtained if the players follow the deterministic strategy denoted by 
(OO^Ol)” - ^, k = 0 , • ■ ■ ,n. Here, for a fixed k the previous notation means that the first 
k players answer aways the output 0 while the next n — k players answer the output 0 if 
they are asked question 0 and they output 1 if they are asked question 1. Interestingly, 
one can show in addition that only two of these strategies (k = 0 and k = n) are relevant 
in average. Indeed, 0 Theorem 4] and |3) Theorem 5] show, respectively, that 

(3-2) p| max{|(T, 7 o)|,|(T, 7 n )|}| = °- 8475 -^+ °( 1 ) ^ 

^ ' n 4 


and 

(3.3) 





for any constant c > 0. 


19 r* 

We will not talk about symmetric correlation Bell inequalities since these should be defined as those 
inequality which are invariant under permutations of the parties without any extra restriction. 
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A simplification in the analysis of the classical value of a symmetric XOR game allows the 
authors in [3j to prove Eq. (13.211 and Eq. (13.31) by using classical probabilistic techniques. 
One can deduce from here that for every e > 0 there exist positive constants 01 (e), 02 (e) 
such that the probability of the classical value w(T) being in [Oi(e)n^ 1 / 4 , 02 (e)n -1 / 4 ] is 
at least 1 — e. This is done by using that the random variable max {| (T, 70 ) |, | (T, q n ) |} 
converges (as n tends to infinity) to the sum of two Gaussian random variables with 
known mean and variance (see [21 Section 5.2] for details). However, note that this does 
not allow to state the previous estimate for fixed constants Oi and O 2 and probability 
1 — o(l). This lack in the concentration of measure seems to be due to the nature of the 
value co(T) rather than to the analysis performed by the authors. 

On the other hand, it was proved in 0 Theorem 1, Theorem 2] that 


(3.4) 


lim P 

n—>• 00 


Ci 


\/In 


n 


1 

n4 


<uo*(T)<C 2 


\/Tnn 

T 

n* 


= 1 for certain constants Ci,C 2 > 0 . 


We see that this result is, in terms of concentration, stronger than the previous one 
for the classical value of the game, since here we can fix the constants C \, C 2 and get 
probability 1 — o(l). To show Eq. (13.411 the authors first used a previous result in [49] . 
[S3] stating that the quantum value of a symmetric XOR game is given by cj*(M) = 


maxi 


1 * 1=1 




Then, computing w*(M) reduces to the maximization 
of the absolute value of a polynomial in one complex variables. On the other hand, the 
previous expression is trivially upper and lower bounded by functions depending on the 
real and imaginary part of the corresponding polynomial, 


max 

ae[0,27t] 


1 

jEH) 

J—0 


Mi 


( n .) COS (ja) 

, max 

\jJ 

aE[0,27r] 




Mi 


3=0 


n 


sin (ja) 


If (Mo, Mi, • • • , M n ) are chosen at random, these expressions reduce to random trigono¬ 
metric polynomials studied in [42] and the problem is actually reduced to study the 
quantities M n (t) = max ae [ 0i27r ] \P n (x,t)\ for P n {x,t ) = Xm=o r m^m{t) cos mx, where 
is the Rademacher system. Although this is the key object of study in [42], the 
results in this work do not fit directly in the problem considered in [ 2 ] and obtaining the 
estimate ([H requires a careful study of the particular case. 


3.2. Sampling quantum n-partite states. A different way of studying Bell violations 
consists of sampling quantum states according to some random procedure and looking 
at the violations they can produce. Note that this is not the same as sampling quantum 
correlations at random as it was considered in Section [TJ In the current situation one 
should regard both the POVMs and the Bell inequalities as free parameters in the 
problem. One can then define a natural measure of how nonlocal a quantum state is 
(see [35]). However, the great freedom in the problem makes it difficult to handle. An 
intimately related problem is to study the probability of a quantum state being nonlocal 
with respect to a particular family of Bell inequalities. This is the context considered 
in [ 22 ], where the authors restricted to the setting of n-partite binary inputs correlation 
Bell inequalities (T Xlj ... ,x n )xi,- ,£„€{(),!}• ^ is well known that in this case the set of 
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correlation Bell inequalities can be replaced by the single nonlinear inequality |50j , 

m 


E 


E 

l<i<n a?iE{0,l} 


n 

3 = 1 


Al + (-1 pA{ 


< 1, 


where Aq,A\ = ±1, 1 < j <n represent the deterministic measurement results for the 
pair of measurements in site j. Hence, one can define the function LV : S 2cfI_1 —)• [R+ on 
the unit sphere of (C d )® n ; that is, the set of pure quantum state \i/)) of n d-dimensional 
systems, by 


LV{\ip)) : = sup 


E E 

l<i<n XiG{0,l} 




4 + Ij-lpAj 


3 = 1 


■0 


Here, the supremum runs over all possible choices of pair of observables A J 0 , A\ per site 
j = 1, • • • ,n. The main result in [22] describes the probabilistic behavior of LV with 
respect to the uniform measure on S 2dn_1 . 


Theorem 3.1. Let n and d be two natural numbers larger than 1, let \if) be a unit vector 
distributed according to the uniform measure in the unit sphere of (<C d )® n and 

denote A v = {\f>) ■ LV\^) > v}. Then, the following inequality holds: 


(3.5) 


P(A) < 2 


n2 n+l d 2 


+ 2 


2 d 2 r 


(v-5-c d 


9 7T^ 


\2, d\r. 
i i) 


Here, 6 is any positive number, v > c^n + 5 and Cd^ n = (|) 2 + 


As we mentioned in the previous section, for binary inputs correlation Bell inequalities 

n 

the optimal quotient between the quantum and the classical value is 2^. Note that Eq. 
(13.51) implies that for large n most pure states do not get even close to this violation. 
For the particular case d = 2 (qubits) one has that Cd t n = 1 f° r every n, so as long as 
v > cn for a suitable constant c, P(A„) -A 0 as n goes to infinity. It was conjectured 
in [39] that this is also true if v is of order y/n log n but this problem seems to remain 
open (see (221 Section I] for details). On the other hand, if d > 3 Eq. (13.51) implies that 
most of the states do not violate any of these inequalities for n large enough. Indeed, 
note that lim n Cd : n = < 1, so we can choose 5 > 0 with 5 + Cd t n < v < 1 such that 

lim n P(A„) = 0 super exponentially. It could be interesting to study this problem when 
one samples on not necessarily pure states. Although the way of sampling in that case 
is not so obvious, there are very interesting works showing different behaviors in that 
picture (see m , m, m and the references therein). 

Theorem 13.II is proved via a concentration-type argument. First, the authors consider 
an e-net for the set Q of possible pairs of observables Ag, A\ per site j = 1, • ■ ■ , n. This 
allows them to fix a finite set Q e of these elements to work with. On the other hand, 
for one such a choice Q € Q e , one considers the corresponding function LV^(\f>)) : 
S 2cfl —i -A P + . It turns out that this function is quite regular in terms of its Lipschitz 

l^Here, we really mean the equations defining the facets of the set C. 











RANDOM CONSTRUCTIONS IN BELL INEQUALITIES 


19 


constant. Then, one uses this regularity in two different ways. First of all, it allows to 
pass from taking the supremum on Q to considering Q e by means of the inequality 

P(A,) := P{|i/>) : sup LV®\ip) > u} < P{|i/’) : sup LV®(\ip)) > v — 5}. 

QeQ Q^Qe 

Secondly, for a fixed Q € Q e one can compute the expectation \E[LV^] and use Levy’s 
Lemma to bound the probability of the points for which the function is far from its 
expectation P(|/ — E[f] | > e) < 2 exp ( — ^ 3 ^ ), where here / : S n —>• P is a real 
function and A is its Lipschitz constant with respect to the euclidean distance. Putting 
the two previous ingredients together allows the authors to obtain Eq. (13.51) by using a 
counting argument. 

3.3. Sampling random measurements. In [30], [43] . [ 35] , the authors considered 
quantum correlations sampled in such a way that the quantum state is fixed and the 
observables are chosen at random. This seems to be a very interesting problem from an 
experimental point of view since it is motivated by the requirement of well calibrated 
devices and a common reference frame between the parties in the standard Bell exper¬ 
iments. Showing that random measurements produce Bell violations implies that the 
previous assumptions can be removed. In [30], [38] the authors considered the n-partite 
2-dimensional GHZ state: | ip) = 2 -1 / 2 JT = 1 1 i)® n and studied the correlations when each 
party j measures with two possible observables Aij, A 2 J chosen at random according 
to the law (RIM): Ajj = mj ■ a, where are uniform and independent vectors on 
the unit sphere of C 3 and a = {(Jx^cTy^Gz) is defined via the Pauli matrices. Numerical 
computations suggest that the corresponding correlations are nonlocal with probability 
tending to one as n goes to infinity. This behavior is also studied when the two mea¬ 
surement per party are chosen so that n\ and are orthogonal but still random (ROM) 
to show that nonlocality is even more likely in this last case. In [33] the authors also 
showed some numerical evidences that in the bipartite case, correlations produced by 
using N random measurements on each party of = 2 " 1 / 2 ]P 2 =1 |i )® 2 lead to a similar 
behavior to the previous one. Since the techniques in these papers are a bit different 
from those considered in the current survey we will not explain them in detail. However, 
we refer the reader to the original works, where many other results can be found as well 
as some very nice explanations about the experimental interpretation of the different 
ways of sampling. 

4. Unbounded violations for bipartite Bell inequalities 

The setting of general Bell inequalities is much richer than the context of correlation 
Bell inequalities. The new parameter given by the number of outputs K leads to study 
the quantity LV n (N,K,d ) analogous to LV n (N,d ) introduced in Section [T71 where n 
is the number of parties. The application of the parallel repetition theorem [311] to the 
magic square game (or any pseudo-telepathy game El) and also the work [29] show 
that LV 2 (N,K,d) cannot be upper bounded by a uniform constant independent of the 
parameters N, K and d (see |26l Introduction] for details). That is, there exist unbounded 
violations of bipartite Bell inequalities. However, the previous unbounded Bell violations 
are far from the upper bounds provided in [26l Theorem 14]. 
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Theorem 4.1. The following upper bound holds: 

LV 2 (IV, K , d) = O (/i), where h = min {IV, K, d j. 

In the paper m the authors used a random construction to show the lower bound 
LV 2 ([ 2 log ^/ 2 ]LlV,IV)=n(\/]V/log 2 IV), improving the previous estimates so far. Since 
this result has been recently improved in two different ways we will omit the details in 
this survey. In the current section we will explain a more recent result which improves 
and simplifies the work in m- Based on a random construction, in [26] the authors 
proved that LV 2 (N,N,N) = fl(\/lV /y/logN). According to Theorem 14.11 the previous 
estimate is only quadratically off from the best upper bounds in all the parameters of 
the problem at the same time. Actually, the most relevant point of this result with 
respect to the previous ones is the decrease in the number of inputs from exponential to 
polynomial. We will also explain in Section B~2l the existence of a bipartite Bell inequality 
providing the lower bound LV2(2 n /N, N, N) = Q(N/ log 2 IV). This gives an essentially 
optimal upper bound as a function of the number of outputs and the dimension of the 
Hilbert space. 


4.1. Unbounded violation with polynomially many parameters. In [26] the au¬ 
thors introduced the following construction: Consider a fixed number of ±1 signs e x a 
with x, a, k = 1, ■ ■ ■ , N. For a constant K one defines the vectors \uf) = ^-^ (1, e x x a . ■ ■ ■ , e x a ) 
in C N+1 for every x,a = 1, • • • ,1V; and 


a) Bell inequality coefficients: (Mx’,y) x ^^ =1 given by 


1 


N 


l\fa,b _ W f k 
x,y jy 2 /_ c x,a*-y,b 

k= 1 

a,b _ 


for x,y,a,b = 1, • • • ,N and Mf" y = 0 otherwise, 

b) POVMs measurements: in My+i as 


K 


l^)<*“| for a = 1,- • • ,1V, 

1-EL xK for a = IV + 1 


for x = 1, • • • ,1V. 

c) States: Let | ip a ) = ELL where (oj)Ei 1 be a decreasing and positive sequence 

verifying ELL a l = L 


The following result can be found in [26, Theorem 2]. 


Theorem 4.2. There exist universal constants C and K such that for every natural 
number N there exists a choice of signs {eLlL, fc=1 verifying that {U“}L=L define 
POVMs measurements, u>(M) < Clog N and 

N,N+1 N+l 

® E h y \y a ) > ai ai • 

x,y-,a,b=l z=2 


(4.1) 
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Moreover, the probability of the elements (choices of signs) verifying Eq. when they 

are chosen independently and uniformly at random tends to 1 exponentially fast as N 
tends to infinity. 

According to the previous theorem, by considering an asymmetric element of the form 
= l/\/2 and ag = l/\/2N one obtains the estimate LV{M) = Vt{y/N /log IV) 0- As 
in the case of tripartite bell inequalities, a probabilistic construction provides a good 
lower bound for LI/^IV, N, N ) which almost matches the known upper bounds. 

The estimate lo{M) = O(loglV) in Theorem 14.21 can be obtained from well known 
probabilistic results. The original proof in [26] was based on Chevet’s inequality but it 
was improved and simplified in SH leading to u(M) = 0(y/log N). This last proof is 
based on standard estimates on gaussian variables. 

The key point to understand Eq. (ED is to consider the families of rank one operators 
E x = ^(KXOlCi for every x, where < = ^=(4,a>"- ,e£ a ). Then, it is not 
difficult to see that if one defines the rank one operator g = |00)(i/’|, where \if) = 
K*)> one obtains 

N 

(4.2) £ K;itr(E° ® E y bV ) = n(y/N). 

x,y;a,b= 1 

The problem here is that the elements (Ef) xa=] do not dehne a family of POVMs (the 
operators are not even positive!) and that the operator g is not a state. Nevertheless, the 
previous family of operators should be understood as a family of non-positive POVMs. 
The interesting point here is that the good properties of these operators assure that by 
replacing the vectors tt“ by u% as in the statement of Theorem 14.21 one obtains a new 
family of operators Ef = |u“)(u“| verifying ^/) a=1 Ef. < 1. The new elements Ef are 
defined by adding more entries to the matrices defining the operators Ef so that they 
are positive and verify the required property on their sum. In particular, Ef is encoded 
in the first column of Ef up to an extra element equal to 1. Then, the reader can guess 
and easily check that if one considers the state \ip) = l/\/2(|00) + IV’)) , the terms |00)(i/;| 
and |i/’)(00| appearing in p = |<^)(<^| have the same effect as in Eq. (14.211 . Indeed, the 
rest of the terms in p do not play any role. Actually, a similar argument can be used for 
a general state |(^ a ) = 1 a i\ti) to obtain the estimate (14.11) (see (20 Section 3] for 

details). 

Finally, the reason to pass from M to M is that the the previous argument gives an 
element (tr(E% <8> E y b p)) x y a b=1 for which Ea=i K < 1 for every x. Hence, by naively 
adding an extra output a = N + 1 per measurement x, the previous modification on the 
element M allows one to complete the family of operators to obtain E x = 1- 

Interestingly, if one plugs the (N + l)-dimensional maximally entangled state \if) in 
Eq. (14.Ill , no large violation is obtained. In fact, it is an open question whether one can 
get large Bell violations on M by using the maximally entangled state, even if we do not 


14 It was shown in [4T] that oj(M) = Of-^/log N) in Theorem 14.21 leading to LV(M) = Q(y/N/ log N). 
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restrict its dimension. In [26] the authors made a modification to the inequality M to 
obtain such an example. More precisely, [26l Theorem 3] shows the existence of a Bell 
inequality M with 2 N inputs and IV +1 outputs par party, and POVMs {E^} x , a acting 
on C^ +1 such that M and {E^} X)a verify the same properties as those in the statement 
of Theorem 14. 2 1 and, in addition, sup{|(M, Q m ax)\] = O (log IV), where this sup runs over 
all quantum probability distributions Qmax constructed with the maximally entangled 
state in any dimension. Previous results in this direction had been obtained in m, m, 
where the authors had proved that certain quantum probability distributions cannot be 
written by using the maximally entangled state. We will not explain [261 Theorem 3] in 
detail since its proof is quite technical and it is based on proving certain estimates on 
completely bounded norms. Furthermore, an improvement of the previous construction 
was made in mi Section 3], where the author gave an explicit Bell inequality with 2 N /N 
inputs and N outputs per party verifying the same properties as M. The proof of this 
result is based on deep techniques from quantum information theory. 

As in the case of Theorem 12.51 one can understand the statement of Theorem 14.21 as 
a result providing large Bell violations with high probability when the Bell inequalities 
are properly sampled. However, while in the case of tripartite correlations the sampling 
procedure was quite artificial, involving in particular a duplication of indices, in the 
current situation the sampling looks quite natural if one ignores the added extra zeros, 
which do not play an important role. It is not clear for the author if some other samplings 
should be considered more (or less) reasonable. 

The lack in the sampling criterium becomes more extreme when one samples quantum 
probability distributions. This has motivated the works in this direction to focus on very 
particular situations. An example of this can be found in [6], where the authors study the 
violation of the CGLMP bipartite Bell inequality M ggmlp with binary inputs and K 
possible outputs per measurement [20]. Actually, even though the inequality M ggmlp 
is fixed, the authors in [6] need an assumption about the best POVMs for that inequality 
to be able to analyze its probabilistic behavior. Once the measurements are fixed, the 
problem of sampling quantum probability distributions is reduced to sampling quantum 
states. The main result in [6] gives the expected value of the violation for the Bell 
inequality M GGM lp w ith respect to the uniform measure on li-dimensional bipartite 
pure states; that is, when the pure states are uniformly sampled from the unit sphere 
of C K . This is done analytically for K = 2 and in the range of large K, while for 
intermediate values it is done numerically. The key point of this study is that the previous 
restrictions allow the author to write the quantum value of the corresponding inequality 
as a function of the Schmidt coefficients of the states (see El Eq. (3)]). Interestingly, 
even for this restricted setting the authors need to use nontrivial techniques from random 
matrix theory to analyze the problem. 


4.2. A comment about some sharp explicit constructions. We will finish this 
work by briefly mentioning two fully explicit constructions introduced in m ■ Since 
this survey is focused on random constructions, we will not go on these examples in 
detail and we refer the reader to the original source. However, the relevance of these 
examples forces us to mention them. The first one, called Hidden matching game , gives 
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a bipartite Bell inequality with 2 N inputs and N outputs per party, which leads to a 
violation of order v/iV/ log N. Although this order is slightly worse than the one in 
Theorem 14.21 and some of the parameters are exponentially higher, this Bell inequality 
presents a very particular form. Indeed, it is introduced via a two-prover one-round 
game G for which uj*(G) = 1. However, we must point out that the authors in [17] 
show that the quotient between the quantum and the classical bias of the game G is 
(3* (G) / (3(G) = Vt(\fN/ log N). From here one can obtain a Bell inequality M (which 
is not a two-prover one-round game anymore and it has different quantum and classi¬ 
cal values from G) such that u*(M)/u(M) = Q(^/N/logN) (see m Section 2] for a 
discussion about these things). The second example provided in |T7] is the KV game , 
which was first introduced by Khot and Visnoi in [33J to show a large integrality gap 
for a SDP relaxation of certain complexity problems. In m the authors carried out a 
careful analysis of the game to show that it provides a Bell inequality with 2 N /N inputs 
and N outputs per party, leading to a Bell violation of order N/ log 2 N. Since this order 
is attained on a quantum probability distribution constructed with the IV-dimensional 
maximally entangled state, it follows that LV2(2 N /N, N, N ) = Q(N/ log 2 N). According 
to Theorem ED the previous result is essentially optimal in both the number of out¬ 
puts and the dimension of the Hilbert space. Finally, this Bell inequality is particularly 
interesting because it is a two-prover one-round game, which implies some additional 
structures (see m for details). As an interesting remark we should point out that the 
fact that the maximally entangled state is the key element to study the KV game is not 
a coincidence. Indeed, it was proved in [26l Theorem 10] that for any Bell inequality 
with positive coefficients (in particular, any two-prover one-round game) the maximally 
entangled state always gives the largest violation up to a logarithmic factor in the di¬ 
mension of the Hilbert space. In this sense, the KV game is very different from the Bell 
inequality given in Theorem 14.21 

It would be interesting to study if one can obtain a similar violation to the one given by 
the KV game by reducing the number of inputs (at least from exponential to polynomial). 
Moreover, there is still a gap of order between the best lower and upper bounds for 
the bipartite Bell violations as function of the number of inputs N. Finally, the quantity 
LVs(N, K,d) seems to be completely unexplored. One can check that some examples 
such as the one introduced in Theorem 14.21 can be generalized to the multipartite setting 
to provide nontrivial lower bounds. Finding an analogous inequality to the KV game in 
the tripartite case seems to be a more challenging problem. 
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